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SUMMARY 


A one-dimensional theoretical method for the prediction of the internal 
performance of a scramjet engine is presented. Using this method, the evalua- 
tion of the effects of vehicle forebody flow parameters and characteristics on 
predicted thrust for the Langley scramjet engine indicate the engine capture 
mass flow, and thus thrust is strongly influenced by changes in vehicle fore- 
body characteristics. Evaluation of the effects of changes in the engine 
internal parameters on thrust indicate that changes in engine flow parameters, 
such as additive drag, fuel equivalence ratio, nozzle length or nozzle exit 
area and combustor pressure integral cause significant changes in internal 
thrust. Theoretical internal performance predictions, in terms of thrust 
coefficient and specific impulse, are presented for free-stream Mach numbers 
of 5, 6, and 7, free-stream dynamic pressure of 23,940 N/m 2 , forebody angles 
of 4.6° to 14.6° and a fuel equivalence ratio of 1.0. 


INTRODUCTION 


Vehicle design considerations in hypersonic flow (reference 1) show that \ 
the reliable prediction of a scramjet engine's performance is an absolute \ 

requirement. This becomes evident upon consideration of the fact that the net ‘ 
thrust is the difference between two large and nearly equal forces, the nozzle 
thrust and the inlet ram drag. Misdirection or location (relative to vehicle 
center of gravity) of the large nozzle thrust vector can result in exorbitant 
vehicle trim drag penalties (reference 1) while Intelligent engine vehicle 
integration to control the thrust vector location and orientation can result 
in significant gains in lift (references 2 and 3). In order to further study 
the problems inherent In hypersonic vehicle engine integration, the need exists 
for engine performance computations covering a range of Mach numbers and angles 
of attack. As any initial hypersonic flight experiments will probably be 
conducted on a small scale experimental vehicle, emphasis should be placed on 
providing hypersonic engine performance results which would be applicable to 
this case. 

To meet this need, the Langley scranjet modular engine (figure 1) discussed 
in references 4 through 8 was chosen fo" analysis of its performance character- 
istics. The modular engine concept lends Itself to integration with the 


ii Mi at 


r 


m i - j tu 4 e • no Inc was designed with modest contraction ratios and 

mid struts are used to enhance the 

inlet flow compression and reduce the combustor length by providing multiple 
injection* - planes? The engine configuration Is 

chart!* has sweat sidewall planar surfaces, has an opening upstream of the cowl 
leading edge through which flow spillage occurs for stortl ^ and noriw 1 i 

Hnn and has the external cowl surface at a small angle (1 1/2°) with the local 
f^\ikeepexterSalwi drag minimal. Vehicle-engine integration Is assumed 
Lch thS l^e portions of the vehicle afterbody are used as part of the 

engine exhaust nozzle. 

toc.K4^ 

SIS SS&S V An ders00 

si deration The cycle process and performance parameters are based on assuming 
one-dimensional fllid flow and real gas eqonibrto the^yn^icpropert^es 
• rnn i..nrtinn with the control volume concept. The thermodynamic property 
Ha t.a^was^ taken film Iefe™« 12 and modified so that the enthalpy base matches 
that of reference 13 The original computerized method has been modified so as 
telcctelf tXIduIi ly and°separatel y for heat (due to surface heat transfer) 
and friction losses of the flow before entering the engine 

^; o ^hr a L^rsfr^^^?I^TS^islI^o^ s l??4clr , o^t-sp?^er h ;^rh t ^hr 

bottom of the Inlet and of Injected fuel temperature, ^ortkinetic energy 
efficiencies are used to account for nonviscous losses of the en 9 ] ne ^apyure 
flow in the forebody, inlet and nozzle portions of the engine cycle analys . 

The method and the modifications are discussed. 

A cycle and performance analysis has been performed for *" J"? 1 J® ' ™^!f c 
of a size believed to be representative of that for an experimental hyp 
vehicle. These computations were made for free-stream Mach number* of 5, 6, 
and 7, for altitudes corresDonding to a free-stream dynamic pressure o 

• 2 ^ 40 f% ( \ 5 ?o 0 TAzW 'SdW It^was C also C tecessary°to 3?" 

estimates of friction and heat losses, nozzle plume drag and size, inlet spill- 
Me IddUlve drag, drag due to lift, end Inlet contraction ratios The 
procedures utilized to obtain these estimates are discussed and the results a 
presented along with the performance results. 


SYMBOLS 



p 

area (m ) 

the area at station "1" that the mass flow VI (no spillage) 
would occupy with the station "l" velocity a8d density 
which includes boundary layer losses (m‘) 
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Mfuel 
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Pr 
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combustor wall pressure Integral parameter defined by 
equation (47) 

engine force coefficient 

engine force coefficient (Ag/q 0 A c ) 

nozzle plume drag coefficient (Dg/q^) 

friction coefficient 

inlet additive drag coefficient (D^/q Q A c ) 

nozzle plume drag (N) 

inlet additive drag (N) 

defined by equation (16) 

defined by equation (15) 

defined as equal to p (N/m) 

friction losses (N) 

acceleration due to gravity (m/sec ) 

enthalpy ( J/Kg) 

heating value of fuel (J/Kg) 

engine height (m) 

engine component heat loss (J/sec) 

heat added to fuel through regenerative engine cooling 

specific Impulse (sec) 

ratio of specific heats 

Mach number 

pressure (N/m 2 ) 

Prandtl number 
dynamic pressure (N/m 2 ) 

Reynolds number based on momentum thickness 
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Subscripts 

A 


entropy (J/®IC) 
temperature (°K) 
increment of temperature (°K) 
velocity (m/sec) 
mass flow (kg/ sec) 

mass averaged nozzle exit flow direction relative to 
vehicle forebody surface 

parameter defined by equation (7) 

combustion efficiency 

shock kinetic energy efficiency (adiabatic) 

component in the flight direction of the engine force 
parallel to vehicle forebody surface (N) 

component in the flight direction of the engine force 
normal to vehicle forebody surface (N) 

absolute pressure force on the capture flow stream tube 
over the inlet spillage flow area and normal to forebody 
surface and defined by equation (3) (N) 

absolute pressure force on the lower surface of the plume 
normal to forebody surface 

parameter defined by equation (9) 

density (kg/m^) 

parameter defined by equation (8) 

engine thrust (N) 

fuel equivalence ratio 

parameter defined by equation (6) 

vehicle forebody angle with free-stream (deg.) 

momentum pA + WV, (N) 


parameters include no boundary layer losses 
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AV 

c 

C£ 

CO 

comb 

f 

fuel 

6 

HS 

1 

inv. 

M 0 = 6 
N 

PE 

S 

t 

W 

0 

1 

2 

3 

4 


Average of parameter over engine component being considered 
cowl 

value of parameters at entrance of engine component being 
considered 

engine component being considered such as vehicle forebody, 
Inlet, combustor, and nozzle 

combustor engine component 

vehicle forebody engine component 

fuel 

guessed value during Iteration for solution 

parameter value as computed for heat sink engine of 
reference 18 

Inlet engine component 

Invlscld value outside boundary layer 

value of parameter at freestream Mach number of 6.0 

nozzle engine component 

Langley scramjet engine module of present paper. Capture 
area Is 0.1672 m2 

Isentroplc expansion of flow at an engine station to 
pressure of the Imnedlate upstream engine station 

stagnation value 

wall value 

free-stream station 

engine Inlet entrance station 

engine Inlet throat station 

engine combustor exit station 

engine nozzle exit station 


Superscripts 

1 prime quantities Include no stock losses for engine component 

In consideration 


ENGINE PERFORMANCE ANALYSIS 


The present method of solution for a scramjet consists of combining a one- 
dimensional fluid flow model and real gas equilibrium thermodynamics with the 
control volume concept. The control volume assumed relative to the engine 
inlet entrance station and the one dimensional fluid flow model chosen are 
illustrated in figure 2. The engine capture mass flow passes from the free 
stream through the vehicle bow shock, along the forebody surface, through the 
inlet (spilling flow through the bottom of the inlet), through the combustor 
and out through the nozzle exit. The changes In the engine flow parameters 
through the engine are computed one dimensionally. Calculations are done 
separately for the vehicle forebody flow field. Inlet, combustor, and nozzle 
and approximately account for the shock losses, friction losses, and heat 
losses along with flow spillage through the bottom of the engine. 

The resultant force (or thrust) in the flight direction exerted on a 
vehicle by the engine is the Integral of the absolute pressure forces In this 
direction on the surfaces washed by the capture flow of the engine. Therefore, 
utilizing the control volume concept and the one-dimensional flow model 
(figure 2) the engine's resultant thrust Is predicted by the computation of the 
engine force A^/cos w in a direction parallel to the vehicle forebody surface 

and the engine force A g /s1n u on the capture flow stream tube In a direction 

normal to the vehicle forebody surface and vectorally resolving both forces to 
the flight direction. The vectoral relationship between the engine forces 
A^/cos w and A g /sin w with resultant thrust In the flight direction Is shown 

in figure 2. The force A^/cos w Is equal to the Integral of the absolute 

pressure forces on the surface washed by the flow and In a direction parallel 
to the forebody surface. The force A^/cos w Is also equal to the nozzle exit 

flow momentum, minus the momentum of the capture flow at the Inlet entrance 
station, the plume drag, and the additive drag sunned in a direction parallel 
to the vehicle forebody. The force A g /s1n w Is equal to the Integral of the 

absolute pressure forces on the surfaces washed (excluding the vehicle forebody 
surface) by the capture flow and In a direction normal to the vehicle forebody 
surface. The force A g /sin w is also equal to the pressure force (Aj) on the 

capture flow stream tube over the spillage area of the inlet, the pressure force 
(Ap) over the plume lower surface area and the nozzle exit flow momentum summed 

In a direction normal to the vehicle forebody surface. The forces A g have to 

be evaluated as they have a component In the drag direction (figure 2) which 
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must be Included In the net thrust determination. The resultant engine thrust* 
In coefficient form and In the vehicle flight direction* Is given by* 




Ua) 


The coefficient relative to the force Is given by* 


C 


A 


A A 

% A c 


Ub) 


and the force A A Is given by. 


Nozzle Exit Inlet Entrance Stream Tube 

Momentum Momentum Pressure Forces 

I 1 

A a * {'(p 4 + P 4 v|)A 4 * Cosine 6 (p, + P 1 V,)A 1 , - , q 0 A c (Cj + C^ 1 } Cos <o 1(c) 

The coefficient relative to the force Ag Is given by. 


qTJT 

^0 c 


1(d) 


and the force Ag Is given by. 


Stream Tube 
Pressure Forces 


{*Aj + 


Nozzle Exit 
Momentum 

^ +, (P 4 A 4 + W 4 V 4 )^s1n B) sin u> 


1(e) 


Another measure of the scramjet engine Internal performance Is given by 
specific Impulse. The specific Impulse Is the thrust In flight direction 
divided by the fuel flow per second and Is given by, 


I 


sp 



( 2 ) 


In order to compute the thrust coefficient or specific Impulse, the engine's 
corresponding one-dimensional Internal flow has to be computed based on assumed or 
computed values for the nozzle exit area A 4 , the additive drag coefficient Cj, the 

plume drag coefficient C^, the plume force Ap normal to the forebody surface, 

the fuel flow Wf Ue ^, and the engine's Inlet mass flow spillage and aerodynamic 
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contraction ratio A^/Ag. The nozzle eartt area* the additive drag coefficient* 

the plume drag coefficient* the plume force normal to the forebody surface, and 
the engine Inlet mass flow spillage and aerodynaml* contraction ratio must be 
determined experimentally or approximated theoretically. The parameters chosen 
as Input to the analysis are listed In Table I. The methods used to obtain 
these Input parameters are discussed In the following sections. 

The calculations for the engine capture flow are made separately and In 
the following order: free stream, vehicle forebody flow field, inlet, combustor, 

and nozzle. In general, the vehicle forebody flow field and inlet calculations 
consist of making constant pressure momentum balance computations separately 
across the vehicle forebody flow field and inlet assuming the pressures to be, 
respectively, p Q and p^. In so doing, the capture flow is penalized as to 

friction and heat losses. Following the constant pressure computation the 
engine capture flow is penalized as to shock losses using the concept of shock 
kinetic energy efficiency. Then in the vehicle forebody flow calculation the 
flow is compressed to pressure "p-j"* the pressure of the inviscid flow down- 
stream of the forebody shock. In the inlet flow calculation the flow Is com- 
pressed to the inlet aerodynamic throat area "Ag". The flow calculation 

across the combustor consists of setting up a control volume on the combustor 
as shown in figure 3, and generating a one-dimensional momentum and heat balance 
for the control volume. The resulting set of equations are solved for the com- 
bustor exit flow conditions. The one-dimensional flow computations for the 
nozzle consist of isentropically expanding the flow from the combustor exit to 
the nozzle exit area A^. The nozzle exit flow conditions thus obtained are 

then penalized as to friction and heat losses in the nozzle. A more detailed 
development of the theoretical one-dimensional flow method utilized is presented 
in the Appendix. 


ENGINE PERFORMANCE 


Prediction of the performance of a scramjet engine requires knowledge of 
several parameters whose values are strongly engine design or vehicle-engine 
integration oriented. In fact, for some engines and vehicles, the determination 
of the values of these parameters can be reliably done only by experimental 
means or at least by a combination of theoretical and experimental means. This 
Is particularly true when dealing with the Langley scramjet module and Its 
integration with a vehicle. For example, schedules (as a function of inlet 
entrance or vehicle forebody inviscid Mach number) of the module spillage flow 

W o " W 1 

— — , additive drag (i^ inlet pressure force Aj normal to the forebody 
0 

surface, and inlet aerodynamic contraction ratio ^/A-j), are determined experi- 
mentally or through a combination of experiment and theory. In the present 
analysis, theoretical approximations are obtained for the nozzle exiv area 
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(A^/A-j ), plume drag (C^), plume force A p normal to forebody surface* and heat 

and friction losses of the vehicle forebody, the Inlet, the combustor and the 
nozzle. The nozzle plume shape, nozzle ex1>^area, nozzle exit flour direction, 
plume force A^, and plane drag are strongly dependent on the vehicle afterbody 

shape (nozzle top surface) In relation to the engine. The methods used to 
generate the additional engine parameters needed to make performance compute- • 

tlons are discussed and the resulting engine parameters presented In figures 4 
through 20. A list of these engine parameters Is presented In Table I. Theo- ^ 

retlcal performance computations, using these parameters, were made for free- 
stream Mach numbers 5, 6, and 7, free-stream dynamic pressure 23,940 N/m z and ^ 
vehicle forebody surface angles of 4.6°, 7.1°, 9.6°, 12.1° and 14.6° relative w 
to the free stream. These predictions are discussed and the results presented 
in figures 21 through 23. The sensitivity of predicted engine force ® 

(parallel to vehicle forebody surface) to the engine parameters of figures 4 
through 20 is also discussed and the results presented graphically in figures 
24 through 29. 


Engine Parameters 


Pressure ratio across vehicle fore body shock and the inviscid Mach number . - 
The pressure ratio p 1 /d q across the vehicle forebody shock, the vehicle 

forebody's inviscid Mach number and the inviscid kinetic energy efficiency 
(n K )f across the shock are assumed to be those computed for conical flow 

corresponding to cone angles equal to the vehicle forebody surface angle "u>" 
relative to the free stream plus 0.4° to account for boundary layer on the 
surface. The resulting theoretical pressure ratios, inviscid Mach numbers and 
inviscid shock kinetic energy efficiencies are presented in figure 4 for free- 
stream Mach numbers 5, 6, and 7. This assumption for the inviscid forebody 
flow is a result of experimental testing of representative vehicle forebody 
shapes as well as comparison of theoretically predicted surface pressures 
(reference 14) with cone surface pressures. 

In let mass flow spillage .- The mass flow captured as a function of the 
mass 71 ow that could be captured without inlet spillage was determined experi- 
mentally. These exDerinental results are presented in figure 5 in terms of 

W - W. 

spillage flow ratio -At for an inlet entrance (or vehicle forebody) 

w o 

inviscid Mach number range (Mj j ny ) of 2.3 to 6.0. In general, the procedure 

used to obtain these results consists of obtaining static and pitot pressure 
surveys of the flow downstream of the struts and integrating the mass flow per 
unit area generated from these surveys over the survey area while correcting 
for local flow direction. A detailed discussion of the method used for the 
Mach number 6.0 data is presented in reference 4 while a more general discussion 
of all the data is presented in reference 6. For the present calculations it 
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; 


was also assumed that 


V-A, 


Is equal to 


Vi 


The cowl area A. Is 


the area obtained by the projection of the Inlet leading edges on a plane nor- 
mal to the forebody flow. This assumption results In the computation of W Q 

(eq. 21 of Appendix) as based on the mass flow captured as computed from 

the one-dimensional flow at station 1 (eq. 20 of Appendix). This assumes the 
one-dimensional flow computed for station 1 extends all across the cowl area. 
The W Q Is therefore smaller than would be obtained if nonviscous flow condi- 
tions were assumed across the spillage flow area (A c - A-j ). 


Inlet aerodynamic contraction ratio .- The Inlet aerodynamic contraction 
ratio was determined experimentally (reference 6) and the results are presented 
in figure 6 for an Inlet entrance (or vehicle forebody) Inviscid Mach number 
range of 2.3 to 6.0. This Mach number schedule Is based on experimental values 
for mass weighted average throat Mach number and total pressure along with the 
assumption of adiabatic flow. 


Inlet normal force and additive drag. - The inlet normal force Aj is the 

absolute pressure force on the capture stream tube over the spillage area and 
is assumed to be approximated by. 




Aj = 2.4 H 2 p 1 (3) 

The parameter H is the engine height. The relation thus assumed for Aj is 

based on the assumption that the lower surface cf the inlet stream tube exposed 
to the external flow and located between the inlet entrance station and the 
cowl leading edge is acted on by the forebody flow field static pressure p, 
over the exposed stream tube surface area. 

The additive drag is determined using a combination of experimental and 
theoretical results. Experimental values for engine capture flow (modified 
values of reference 6 between (Mj) Inv values of 2.3 to 6.0) are used to deter- 
mine approximately where the dividing stream line (between capture and spillage 
flow) is located at the inlet entrance station (station "1") to the cowl lip. 
This dividing streamline is determined by matching pressure and flow direc- 
tion along the stream tube surface dividing the Internal and external flow. 

This method approximates the down turning of the flow and the pressure in each 
of the shock bays open to the external flow. This procedure is applied through 
the inlet along the dividing stream line to a point just foreward of the cowl 
lip. The strong shock located in the vicinity of the cowl lip and detached 
shock waves (for the lower Mach numbers) make It necessary to assume a relative- 
ly arbitrary fairing from the dividing stream line, already defined, to the cowl 
lip. The absolute static pressures along this stream line are then approximated 
by the experimentally measured inlet sidewall and strut static pressures. The 
additive drag as approximated by this procedure is presented in figure 7 in 
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terms of additive drag coefficient (C t » D./q. A„)^or free- stream Mach numbers 
of 5, 6, and 7. 1 A 0 c 

♦ 

Nozzle plane drag, normal force A p , exit area and exit flow direction . - 

The Internal nozzle flow and external plume are computed theoretically using a 
technique known as floating shock fitting developed by Morettl In reference 15 
and discussed further by Salas In reference 16. The computerized version of the 
method used In the present nozzle-plume computations was developed by Salas for 
two-dimensional and axlsymnetrlc flows. The two-dimensional version is used In 
the present computations. The Internal nozzle flow properties, external plume 
shape and flow properties including local flow directions, and pressure distri- 
bution along the lower surface of the external plume are generated using the 
nozzle vehicle afterbody relationship presented in figure 8. The nozzle-plume 
computations are for a perfect gas with a specific heat ratio of 1.25 and 
initial conditions at the nozzle entrance are obtained using the method of the 
present paper to generate one-dimensional flow parameters at the nozzle entrance. 
The plume drag is then calculated by integrating the absolute pressure forces 
on the lower surface of the plume in a direction parallel to the vehicle fore- 
body surface. The results of these computations in terms of plume drag coeffi- 
cient (C E ) are presented in figure 9 in terms of D E /q Q A c . The nozzle force 

A p normal to the vehicle forebody surface is obtained by integrating the 

absolute pressure forces on the lower surface of the plume in a direction nor- 
mal to the forebody surface. The results of these computations are presented 
in figure 10 in terms of A p /q Q A c . The nozzle exit area is obtained from the 

plume shape at the trailing edge of the vehicle afterbody surface and in a plane 
normal to the average nozzle flow direction. The results of these computations 
in terms of area ratio (A^/A-j) are presented in figures 11(a) and 11(b). The 

average nozzle exit flow direction (6) is obtained by mass weighting the local 
flow directions across the nozzle exit plane and the results of these computa- 
tions are presented in figure 12. « 

Friction and heat losses of engine capture flow. - The friction losses are 
computed assuming the form of Reynolds analogy put forth in equation (2) of 
the Appendix. Therefore, before friction losses can be computed, an approxima- 
tion of the heat loss must be obtained. 

Heat-transfer computations were made for the vehicle forebody region ahead 
of the engine and the results are presented in figure 13. These theoretical 
heat-transfer values for the forebody were computed using the method of refer- 
ence 17, assuming flat plate flow and transition beginning at R 0 = 1000. The 

assumption of flat pla^.e flow for these computations was made as it Is believed 
that this assumption will produce a thicker boundary layer than that for a 
corresponding conical flow assumption and therefore would be the pessimistic 
case. 

Heat-transfer computations for the internal surfaces of the engine are 
reported for a smaller scale heat sink version in reference 18. These 
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computations were made for M Q • 7.0, Mj • 6.0, and $ » 1.0 and produced 
* the following heat loss values for the Inlet, combustor, and nozzle. 

(AHj)jjS « 3.26 x 10* Joules 

(AH comb ) HS - I. 6 ® 7 x 10 5 Joules (4) 

(A^)^ » 5.60 X 10 4 Joules 

The above heat-transfer values obtained for the heat sink engine (subscript HS) 
of reference 18 are scaled up to the present engine (subscript PE) for a free- 
stream Mach number of 6.0 using the following approximating relationship 
(modified from reference 17) between the engine component heat-transfers of the 
two engines. 




* 


^CO^PE 



y 

( Vpe 

q HS 


A 

CE 

( Vhs 

Ve. 


^VhsV 


m 




(h t,o " h w*PE 

( h t,0 - Mhs 



(5) 


This relationship for scaling of the heat load assumes the friction coefficient 
is inversely proportional to (q rt A_)V7 instead of the more correct expression 

O C n yp 1/7 • 

of being inversely proportional to (q A ' ) ' . In the scaling of the heat 

U L 

transfers to very large engines it becomes Important that the more correct 
expression be used. The engine component (subscript CO) heat-transfer results 
obtained using this computational procedure are presented in figure 14 for 
M q = 6.0; trie nozzle heat transfers presented do not •Include the heat transfer 

to the vehicle afterbody portion of the nozzle. Results were obtained for the 
free-stream Mach numbers 5 and 7 using the theoretical curves for engine cooling 
presented In figure 23 of reference 7 to generate values for the ratio of engine 
component heat-transfer to the engine component heat-transfer at M Q = 6.0. 

The results of these computations are presented in figure 15 In terms of 
AH/AH^ 6 g for a free-stream Mach number range from 4.0 to 8.0. The engine 

component heat-transfers calculated (for M = 5, 7) using the curves of fig- 
ures 14 and 15 are presented in figure 16; again the nozzle value does not 
Include the heat-transfers for the vehicle afterbody portion of the nozzle. 
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Remaining engine parameters .- Values for the remaining parameters needed 
in order to make engine performance computations are assumed keeping in mind 
any reasonable physical limits or known vehicle restrictions. These include 
fuel equivalence ratio (4>), combustion efficiency (n com j,)» enthalpy level of 

the hydrogen fuel, and combustor wall pressure Integral factor (given in refer- 
ence 19 and discussed in the Appendix). 

For an operational type hypersonic vehicle it would be desirable to 
achieve cruise at a fuel equivalence ratio less than 1.0 so as to have acceler- 
ation capabilities up to b = 1.0. Engine size is a problem for an experimental 
vehicle ar.d it becomes more important to keep engine size down as much as 
possible by designing tor cruise at 4 - 1.0 and accelerate using fuel equiv- 
alence ratios greater than 1.0. Therefore, the major portion of the engine 
performance computations of the present paper were made for a fuel equivalence 
ratio of 1.0 (ana fl = 5, 6, and 7) with representative engine performance 

computations for fuel equivalence ratios of 9.5 and 1.5 (and M = 6). 

The value assumed "or the combustion efficiency (n.Q^) is 0.95 which is 

believed to be a realistic assumption for fuel equivalence ratios of 1.0. For 
fuel equivalence ratios greater than or less than 1.0, the combustion efficiency 
increases out to what decree it is not known. Therefore, the arbitrary assump- 
tion of 1.0 for the combustion efficiency was made for the engine performance 
computations at fuel equivalence ratios of 0.5 and 1.5. 

The cntha'py 1 eve’ of the fuel is assumed to be 1.415 x 10^ Ooules/kg 
which corresponds to a fuel temperature of 111. 1°K; the enthalpy base is that 
of reference 11. The 'inlet, combustor, and enclosed portion of the nozzle are 
assumed to be >eoenerati vely cooled and therefore the heat removed by cooling 
these engine components is added to the fuel. 

Engi ne pa r ame t ers for 0 f 1.0 .- The forebody flow field and inlet flow 

field are theoretically unaffected 'by fuel equivalence ratio and therefore 
values of pert inen* - performance parameters for the forebody and inlet discussed 
for o - t.r would be applicable to other fuel equivalence ratios. Combustor 
and nozzie heat -transfers, and to some extent friction losses, would be affected 
ar.d therefore wan’d nave to be altered to account for a fuel equivalence ratio 
different from 1.0. Nozzle exit area and exit flow direction would also be 
affected along with nozzle plume drag. The effect cf fuel equivalence ratio 
on these parameters was computed for the heat sink engine of reference 18. The 
results of those c a’: Durations are presented in figures 17 through 20, non- 
dimensional i znd wHh resoect to the value of the parameter being considered as 
computer for 1.0. Check computations for o - 1.5, M = 7.0 and 

cj - 9.6° were me hr for the same oarameters for she scramjet module of the 
present cape-* • for -l! except heat-transfer) , and these res: ts are also pre- 
sented in Mgures 1 '> through -1C. 8a sea on the satisfactory comparison shown 
on those figures, it is assumed that at least for -t = 1.5 realistic engine 
performance cornu u u at ion-' (possibly with the exception of heat-transfer) can be 
made using the curves o r figures 17 through 20 and the value of the correspond- 
ing parameter * -* i.0. 
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Predicted Thrust and Specific Impulse 
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presented^n 0 ^ igures^l ^hrough^ Z^i^termis^f Langley scran >Jet engine Is 

IK T e Aa and "" th ^ 

WiSrS^g ss^'sss inr 0 ? 

putations the engine capture flow*? d . en 9j ne ?* In the forebody corn- 

transfer losses ft thf free-stre?m nrS«,^l Z ?J re1 J* f r e to Action and heat 

would occur at the forebody pressure !?!!, (Appendix) . Physically these losses 

gain. A check into the effect of imDosina°tho f« SU k^i ln « a - d ^ erent entropy 

at the forebody pressure revealed anYcrease in Dredi?tlH 1Cti0n and heat 1osses 
1.5 percent to 7 DPrreni- 4 increase in predicted mass capture from 

to 10 Percen^for thfranorof ?nXd aSe \ n thrust coefficient ofY 9 percent 
putations for the cowl drag indicate that Ye Yw? S H dere< - ’ Preliminary com- 
to 6 percent of the predicted thScJ.Jl! th ! C0W I drag 1s less than 4 Percent 
dimensional effects and end effects ° f magnitude of the three- 

heat transfer of the vehicle afteYdYn^ . know "- u The skin-friction drag and 
the aircraft. It wfs dp^dfd ™ zzle a, ' e charged off to 

transfer of the vehicle'afterbSdJ £-I?2 of A I^ l ?"f ri ? tl0 " drag and heat 
mance computations due to the factY\ : *'nn£ f h vehicle in the engine perfor- 
body portion of the nozzle Is to L con C ^ ^ 3nd ? f what this aft ^“ 
that the heat and friction losSs nf i — ‘ K/Wever » if n ts assumed 
nozzle are of the samf o?der of ! V ? h,ClG ^‘ t - boc V Portion of the 

nal portion of the nozzle! theJ fo^the r U "+* are ^ aS that of the I"ter- 
friction and heat losse would i^easl E ° f flgure 8 ’ nozzle flow 

en 9 me force A, is 

23 ,nd,cate that wjth 

(reference 19) a thermal choice will mcw fS^T ^“T- , 1B n* 9ra l pa C ameter 

body surface angle between 9 6° and 12 ,« ° . ! * " 3 " d a f0re ' 

thermal choking of the combustnrYL]^] j . It: n,ust be noted that the 

of assumptions made in reference 19 to obtai^th®” 6 ^ mensionall y Is a result 
choking limit should not be ISSroreted « l nS* pr f ssure in tegral. This 
real limit. The component 

Ms® Mr «**• « ~ 

an increase in vehicle forPhnHt/^ ' n i fl Tu° wn ^ fl P ur e 23 decreases with 
thrust and specific impuls^dVate fr^ e tL° b f erVat '° nS 0n P redicted engine 
and the presently assumed vehirlp pnni ^ view point of the Langley engine 
be designed to ^ 
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Sensitivity of Predicted Engine Force Parallel to Forebody 
Surface to Predicted and Assumed Engine Parameters 


The performance predictions of a scramjet engine are Influenced to various 
degrees by the predicted or assumed values for the engine parameters discussed 
In previous sections. The relative effects on predicted performance (through 
predicted engine force) of these engine parameters are discussed In terms of 
vehicle forebody parameters, engine Inlet parameters, combustor parameters, 
and nozzle-plume parameters. 

Vehicle forebody parameters .- The effect on predicted engine force A^ of 
increasing the forebody surface angle is demonstrated in figure 24 for M Q = 6.0. 
The force A A and mass flow captured are nondimensionalized by the predicted 
force A ^ and mass flow captured for vehicle forebody angle of 4.6° and pre- 
sented in figure 24 as a function of forebody angle w. From the theoretical 
results presented in figure 24 for M = 6.0, it is concluded that an increase 

in the forebody angle (or vehicle angle of attack) results in a increase in the 
predicted mass flow captured by the engine and a like increase in predicted 
force A ft . 

The effects of forebody boundary- layer losses as well as the engine's 
internal boundary layer losses on engine thrust are introduced indirectly 
through the heat-transfer predicted to pass from the engine capture flow through 
the vehicle's forebody and engine's internal surfaces. The friction losses 
are then approximated using the form of Reynolds analogy presented by equa- 
tion (2) of the Appendix. In order to evaluate separately and concurrently the 
sensitivity of predicted force A^ to assumed values for forebody friction 

and heat losses of the capture flow, three sets of computations were made for 
M q * 6.0 and the results are presented in figures 25 through 27. In fig- 
ures 25 through 27 the force A^, mass flow captured, heat losses, and friction 

losses are nondimensionalized with respect to the corresponding parameters 
predicted for the present engine at M Q = 6.0. 

The curves of figure 25 were predicted with friction losses held constant 
while allowing heat losses to vary from those predicted for M Q = 6.0. The 

effects on predicted force A ft of increasing the engine's internal heat losses 

(with the forebody heat loss held constant) are shown in figure 25(a) to be 

quite small; less than 1 percent for 50 percent increase in engine internal 

heat losses. Of particular note is that the capture mass flow remains constant. 

The forebody, inlet, combustor and nozzle flow heat losses assumed for the 

generation of the curves of figure 25(a) are related in the following manner to 

the = 6.0 values of figure 1: 
o 
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AHf • (AHfJpjr 

AHj - xtAHjJpg 

^comb * X^comb^PE 

^ - x(AHn>PE 

The effects on engine force A A and mass flow captured of Increasing the fore 

body heat transfer (possibly by decreasing the forebody surface temperature) 
are illustrated in figures 25(b) and 25(c). The forebody, inlet, combustor, 
and nozzle flow heat losses assumed for the generation of the curves of fig- 
ures 25(b) and 25(c) are related in the following manner to the values for 
M I ■ 6.0 of figure 1: 


AHf * C(AHf)p£ 

AHj = cUHj) pE 
^comb " ^ AH comb^PE 

“n ■ sfaVfE 


(7) 


An Increase In forebody heat-transfer Is shown to result in an Increase In 
engine force A A and mass flow captured. A 23 percent Increase In forebody 

heat-transfer (change In forebody surface temperature from 666. 7°K to 111.1 °K) 
results In a 4.7 percent to 6.6 percent Increase In engine force A A and a 

3.8 percent to 5.2 percent Increase In mass flow captured. The predicted 
Increase In engine mass flow coincides with the fact that cooling the boundary 
layer flow reduces the displacement thickness and therefore the mass flow 
difference of the boundary layer due to friction losses. 

The curves of figure 26 were predicted with heat losses held constant 
while allowing friction losses to vary from those predicted at M Q = 6.0. The 

effects on engine force A A of varying the engine's Internal friction losses 

(with the forebody friction losses held constant) are Indicated In figure 26(a) 
to be small. A 50 percent increase or decrease in engine Internal friction 
losses results, respectively, In a 2.25 percent to 3.35 percent decrease or a 
2.45 percent to 3.52 percent Increase In the engine's predicted force A A 

while capture mass flow remains constant. The forebody. Inlet, combustor and 
nozzle flow friction losses used to generate the curves of figure 26(a) are 
related In the following manner to the values for M * 6.0: 


i 
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<AF f ) *■ 

(AF X ) « 

^COHlb 

(AF n ) ■ 

The effects of varying the vehicle forebody friction losses on predicted values 
of engine force A ft and mass flow captured are Illustrated In figures 26(b) 

and 26(c). The forebody, inlet, combustor, and nozzle flow friction losses 
assumed in the generation of the curves of figures 25(b) and 25(c) are related 
in the following manner to the values for M Q = 6.0: 

AF f = £(AF f ) pE 

AFj = 5(AFj) pE 

(9) 

^comb ~ ^'^conrtrPE 

af n = C(af n ) pe 

A 50 percent increase in forebody friction loss results in a 19.5 percent to 
22.7 percent decrease in the engine force A^ with a corresponding 17.2 per- 
cent to 21.5 percent decrease In engine mass flow. The engine capture mass 
flow without any boundary layer would be between 13 percent to 18 percent 
greater than the present case. Therefore, there is a limit to the amount the 
friction can be reduced relative to the heat loss and still be physically 
realistic. Therefore, no conclusions will be drawn from the reduction in 
friction portion of the curves of figures 26(b) and 26(c). 

The curves of figure 27 were predicted assuming engine heat losses and 
friction losses vary by the same factor from those predicted at M Q = fi.O. 

This is the more realistic relationship to assume relative to possible errors 
In predicted values. Equations (7) and (9) present the manner the forebody, inlet, 
combustor, and nozzle flow heat and friction losses (assumed to generate 
figure 27) are related to the corresponding values for M Q = 6.0. The varia- 
tions of engine mass flow are presented In figure 27(a) and the variations of 
predicted values for engine force are presented in figure 27(b). The 

curves of figure 27 indicate that the beneficial effects of increased forebody 
heat transfer can be realized only by not Increasing the friction losses In the 
same proportion, otherwise they will be more than offset by the corresponding 


i 


Increase in forebody friction losses. 


°^l)pE 

) » a ^comb^PE 
o(tf M ) PE 
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The theoretical results presented In figures 24 through 27 were obtained 
by varying simultaneously and separately various vehicle forebody (forebody 
angle, friction and heat losses) and engine Internal (friction and teat losses) 
flow characteristics. From these theoretical results It Is concluded that the 
engine's predicted capture mass flow and therefore the engine's predicted force 
A a are strongly Influenced by any changes In forebody flow characteristics. 

Therefore, It Is confirmed that knowing the quantity and physical character- 
istics of the mass flow captured by the engine Is of primary Importance for 
predicting engine performance. The theoretical results of figures 24 through 
27 aTsc suggest some possible forebody design considerations which would Improve 
engine capture mass flow as well as thrust. For example, a conical forebody 
and forebody surface cooling (forebody surface washed by engine capture flow) 
would decrease the boundary layer thickness at the engine Inlet face as well 
as decrease the boundary layer mass flow deficiency. Both of these forebody 
characteristics would increase engine capture flows and therefore thrust. 

Engine inlet parameters .- There are five engine inlet parameters for which 
engine performance sensitivities are known or are evaluated and the results 
presented. Two of these, the friction and heat losses, were evaluated and the 
results discussed along with the forebody flow friction and heat loss results. 

A third is the mass flow spilled. Previous results already presented (fig- 
ures 24 through 27) indicate that engine force values are approximately pro- 

portional to engine mass flow captured. The three remaining engine inlet 
parameters are the additive drag, inlet aerodynamic contraction ratio and inlet 
inviscid kinetic energy efficiency. The additive drag, and spillage mass flow 
could possibly be a function of inlet contraction ratio but at the present time 
the relationship is not known. The lowest value computed for the inlet inviscid 
kinetic energy efficiency is .994 (for (M-,). * 6) and Is assumed for all 

calculations. 1 inv * 

The additive drag given by cos o> nondlmenslonallzed by the predicted 

engine force A^ Is presented for M Q * 6.0 In figure 28 (plume drag Is also 

presented) as a function of forebody angle («). The additive drag distribution 
of figure 28 was generated In the manner discussed In a previous section. At 
forebody angles of about 14° a 50 percent change In the predicted additive drag 
can result in a 7.9 percent change (in the opposite direction) In predicted 
engine force A^. This sensitivity of predicted force A^ to the additive 

drag indicates the importance of an accurate prediction of additive drag. 

The effects on engine force A^ of varying the inlet aerodynamic contrac- 
tion ratio from those determined experimentally is illustrated for M = 6.0 
in figure 29. The engine force A ^ predicted by varying the Inlet aerodynamic 

contraction ratio is presented as a function of the aerodynamic contraction 
ratio Aj/A£ with both predicted force A^ and aerodynamic contraction ratio 

divided by the corresponding value assumed for M Q » 6.0. A 20 percent change 
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In Inlet aerodynamic contraction ratio Is shown In figure 29 to result In a 
3 percent to 4 percent change In predicted engine force A^. 

Engine combustor parameters. - There are seven engine combustor parameters 
for which predicted engine performance sensitivities are known or are evaluated 
and the results presented and discussed. These combustor parameters are the 
heat and friction losses, the combustor wall pressure Integral factor, size of 
the combustor exit area. Injected fuel temperature, combustion chemical eff1> 
clency and the fuel equivalence ratio. 

Predicted engine force A^ sensitivities to friction and heat losses are 

discussed in a previous section. The combustor wall pressure integral factor 
is the integral of the actual combustor wall pressure area distribution ratioed 
to a wall pressure area value corresponding to the assumption that the mean 
combustor pressure is an arithmetic average of the combustor entrance and exit 
values. If the arithmetic average of entrance and exit combustor pressures are 
used Instead of the integral, possible errors In predicted engine force A^ 

values at M Q * 8.0 are shown in reference 19 to be as large as 20 percent. 

Therefore, a correlation for the combustor wall pressure Integral factor 
developed in reference 19 is utilized in the present computations. The size 
of the combustor exit area for the Langley scramjet engine is fixed to be 
64 percent of the cowl area and therefore the sensitivity of predicted engine 
force A A values to the size of the combustor exit area is not considered. 

Theoretical computations of engine force revealed a 5.6 percent Increase in 
A a with an Increase in Injected fuel temperature ( non-regenera tive) from 300°K 

to 1000°K; the feasibility of fuel temperatures of 1000°K is questionable. 
Theoretical computations of engine force A A as a function of combustion effi- 
ciency were made for the heat sink engine of reference 18. These theoretical 
computations which were for M 0 = 7 and $ * 1.0 revealed a decrease of com- 
bustion efficiency from 1.0 to 0.85 resulted in a 14.1 percent decrease in 
engine force A^. 

The effects of the fuel equivalence ratio on engine force A^ are Illus- 
trated In figure 30 for M Q * 6.0. The ratio of predicted engine force A A 

to that predicted for <p * 1 . 0 Is presented as a function of fuel equivalence 
ratio <|>. Also presented In figure 30 Is the corresponding engine force A^ 

distribution computed for the heat sink engine of reference 18 but for M « 7, 

q 0 = 19,496 N/m and combustion efficiency n con) ^ * 0.95. The combustion 

efficiencies for the Langley scramjet engine computations are assumed to be 
0.95 for <j> * 1.0 and 1.0 for * 0.5 and 1.5. The conclusion can be drawn 
from figure 30 that predicted engine force A A values are a strong function 

of the fuel equivalence ratio. For fuel equivalence ratios between 1.0 and 0.5 
the decrease In engine force A^ Is predicted to be 41 percent to 45.5 percent 

depending on the vehicle forebody angle. For fuel equivalence ratios between 
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1.0 and 1.5 the Increase In predicted engine force is 18.2 percent to 

19.5 percent depending on the vehicle forebody angle. The conclusion can be 
drawn from the theoretical curves of figure 30 that the Increase In engine force 
a A a PP #irs off above 4 ■ 1.5 which Is consistent with other 

calculations. 

Engine nozzle parameters. - There are five engine nozzle parameters for 
which engine performance sensitivities are known or are to be evaluated and the 
results presented. These nozzle parameters are the nozzle flow heat and fric- 
tion losses, the average nozzle exit flow angle relative to the vehicle fore- 
body, the nozzle plume drag, and the size of the nozzle exit area A. relative 
to the inlet capture area A^ . 4 

Predicted engine force sensitivities to nozzle flow heat and friction 

losses are discussed In a previous section. The average nozzle exit flow angle 
(B) relative to the vehicle forebody Is shown to be between 0° and -8° (fig- 
ure 12) which Influences engine force prediction through multiplication of 

cos 8 times (smallest value of cosine B Is greater than 0.99) the nozzle 
exit momentum. As the engine thrust Is a small difference between two large 
numbers, the nozzle exit flow momentum and the inlet entrance flow momentum, 
this correction on the nozzle exit momentum can result in significant percentage 
changes In predicted engine force A^. For example, computations for the heat 

sink engine (of reference 18) for M Q = 7.0 and <(> = 1.0 for 6=0° and -6° 

revealed 3.8 percent reduction In predicted engine force A^. The nozzle plume 

drag given by cos u nondlmenslonalized by the predicted engine force A^ 

Is presented for M Q = 6.0 In figure 28 as a function of the forebody flow 

angle (w). The plume drag distribution of figure 28 was generated In the manner 
discussed In a previous section. A 50 percent change In the predicted plume 
drag can result In a 3 percent change In predicted engine force A^. 

The effects of reducing the size of the nozzle exit area on predicted 
engine force A A are illustrated for M Q = 6.0 In figure 31. This reduction 

in size of the nozzle exit area was affected by the reduction of the length of 
the vehicle afterbody portion of the nozzle while keeping the same vehicle 
afterbody surface angle. For H Q * 6.0 and the present engine-vehicle config- 
uration, a 49.5 percent reduction In nozzle length results In a 39 percent 
reduction in size of the nozzle exit area and a 13 percent to 15.2 percent reduc- 
tion in predicted engine force A^. 


CONCLUSIONS 


A one-dimensional theoretical method Is presented for the prediction of a 
scramjet engine's Internal performance. This method of solution consists of 
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combining a one-dimensional flow model and real gas equilibrium thermodynamics 
with the control volume concept. Using this one-dimensional theoretical method, 
the evaluation of the effects of vehicle forebody flow parameters and character- 
istics on predicted thrust values for the Langley scramjet engine Indicate that 
the engine's predicted capture mass flow, and thus the predicted thrust, Is 
strongly Influenced. It Is confirmed that knowledge of the quantity and 
physical characteristics of the mass flow captured by the engine Is shown to be 
of primary Importance for predicting engine performance. Also some vehicle 
forebody design considerations such as high surface cooling and shaping the 
forebody are suggested as they would decrease the boundary layer mass defect and 
thus Increase engine mass flow. The theoretical evaluation of the effects of 
any changes in the engine's Internal parameters on predicted thrust values 
indicate changes in some Internal parameters such as additive drag, fuel equiva- 
lence ratio, nozzle length or nozzle exit area and combustor pressure integral 
cause significant changes in values of predicted internal thrust. Changes in 
other internal engine parameters such as heat and friction losses and plume 
drag result in less significant changes in the predicted values of thrust. 

Theoretical Internal engine performance predictions in terms of thrust 
coefficient and specific Impulse are presented for free-stream Mach numbers of 
5, 6, and 7, free-stream dynamic pressure of 23,940 N/ms forebody surface angles 
from 4.6° through 14.6° and a fuel equivalence ratio of 1.0. These engine 
performance predictions indicate combustor thermal choking at free-stream Mach 
number of 5.0 and a vehicle forebody ar.jle between 9.6° and 12.1°; thir thermal 
choking as predicted one dimensionally is a function of the combustor exit size 
and wall pressure integral assumed. From the viewpoint of the predicted inter- 
nal thrust and specific impulse, the vehicle should be designed to cruise with 
vehicle forebody angles of less than 10°- 
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APPENDIX 


Performance Analysis and Solution Technique 


The method of calculating engine performance consists of combining a one- 
dimensional fluid flow model with real gas equilibrium thermodynamic properties 
and the control volume concept Illustrated In figure 2. The capture mass flow 
passes from the free stream, through the vehicle bow shock, along the forebody 
surface, through the inlet (spilling flow through the bottom of the inlet), 
through the combustor and out through the nozzle exit. The changes in the flow 
properties through the engine are computed by means of the one-dimensional 
conservation equations for mass, momentum, and energy. Calculations for the 
capture flow are made separately for the vehicle forebody flow field, inlet, 
combustor, and nozzle. These one-dimensional computations account approximately 
for the shock losses, friction losses, and heat losses that occur in ea«*h com- 
ponent of the engine. Details of the analysis of each component flow process 
and the methods used to obtain a solution are presented in the following 
sections. 


Free-Stream Flow 


The free-stream conditions ahead of the forebody shock (station "o") are 
computed based on the 1962 Standard Atmosphere Tables (reference 20). In order 
to accomplish the computations of the free-stream conditions using reference 20, 
values for the flight altitude and Mach number are required. Initial guesses 
for free-stream capture area (A Q ) and total momentum of the capture flow 

(p Q A 0 + W Q V 0 ) are then computed from these conditions using an initial guess 

for the capture mass flow (W Q ). These free-stream flow quantities are modified 

In the portion of the forebody flow field computations in which the actual 
capture mass flow is computed. 


Forebody Flow 


The change in flow properties from the free stream across the forebody 
flow field to a station on the surface just in front of the Inlet (station 1) 
is computed using a form of the one-dimensional momentum and energy equations. 
These equations account for forebody flow friction and heat losses as well as 
inlet mass flow spillage. A method is also Included which approximately 
accounts for forebody shock losses. The pressure ratio P-j/P 0 across the bow 

shock, the inviscid kinetic energy efficiency (n^)f across the bow shock, inlet 


mass flow spillage to free-stream capture flow excluding spillage 


(W 0 - W, ) 
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are specified. The forebody Invlscld kinetic energy efficiency (%)« Is given 
by, 

y2 

<«> 

1 » S 

The subscript "S" means these parameters correspond to the expansion of station 
"1" flow Isentroplcally to the pressure at the preceding fi<m||tat1on (which In 
this case is station 0). The prime means no forebody shock wises are included 
In this quantity. The forebody friction losses are computed using the form of 
Reynolds analogy as given by. 
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The momentum relationship assumed for the one-dimensional forebody flow 
Is given by. 


or. 
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The momentum relationship of equations (3) and (4) include no forebody shock 
losses. For s equal p Q , 

F < p o’ p l, s' * p o 


(A4) 


(A5) 
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and, combining equations (4) and (S) gives. 


V S-V <Sf) f ^Avlf 


Utilizing the form of Reynolds analogy as put forth In equation (2) gives. 


,p 7^AV f ( i^> f < Cf: A») f 


Pr T t^y), 


combining equations (6) and (7) gives. 


V 1,S' V 0- J 


pr 7 AH f (V Ay ) f 

^,o*Vf )H l 


T r<»1ng the forebody shock losses on V, - through use of the forebody 
shock kinetic energy efficiency gives, * 5 

1 

V l, S* V 1,S ( "K>f (M 

The energy relationship assumed for the one-dimensional forebody flow Is given 
by. 


’<> + -r ’ V s'nH 1 


V 1 S 4H f 

+ J t^ + Ti7 


V 2 V? , AH- 

h i, s * h 0 + 1 1 wr 


with h^ j and ^ known (p-j g * p Q ) an Iteration on the temperature Is 

performed. The first guess for Tj g Is that It Is equal to T Q and 
succeeding guesses are computed using. 
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(A12) 


Thermodynamic property data are then used to find the values of h" for 
iTj s - AT) and h + for (Tj s + AT), respectively, and the Iteration Is con- 
tinued until succeeding guesses of T-j s differ by less than one degree. The 

thermodynamic property data that are used In the present calculatlon^conslst 
of values for enthalpy, entropy, specific heat, and molecular weight“f each 
species Included In the equilibrium calculations. The data are taken from 
reference 12 and the enthalpy base is adjusted to match that of reference 13. 

The station 1 flow conditions are then determined using the Isentroplc 
flow relations and the station 1 Isentropically expanded flow conditions de- 
fined in equations (1) through (12). The isentropic flow relationship between 
station 1 areas A-j and A-j is given by. 



(A13) 


where A-j j is the area at station 1 that mass flow W Q (no spillage) would 

occupy with station 1 velocity and density (Including boundary layer losses). 

The aret' A-] j would be equal to the cowl area if there were no forebody fric- 
tion and heat losses and no Inlet spillage. In order to find station 1 flow 
conditions, an isentropic compression process from the station 1 isentropically 
expanded flow (subscript 1, S) to station 1 flow Is computed using equation (13) 
in conjunction with a double iteration process. This double Iteration process 
consists of assuming a capture mass flow W-j and Iterating on Aq/A-j^ j until 

a value of A-j/A^ s Is obtained which when substituted Into the Isentropic 
relations produces a value of P^/P 0 within 0.001 of the value specified 
initially. The isentropic relations and the method of solution are as follows: 
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Isentroplc delations 
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A trial value of $ Is assumed Initially and successive values of 
are computed until the new F-j changes by less than 0.001. Initial guesses 
for the values of and T-j are obtained from: 


Pi 
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The thermodynamic property data are then used along with (p-j, + AT) and 

(P-j* T-j) to obtain entropy values of sj and Sy respectively. The value of 
the temperature Is improved by the equation 





(A19) 


Where Sj and sj are the entropies for the states corresponding to (p^ T-j) 
and (p.j, T.j + AT), respectively. Successive values of T-j are computed until 
the new value of T-j computed differs from the old value by one degree. The 


i 




26 


value of the capture mass flow w 1 Is then computed from 

"l ’ "l V 1 A c (’ - (MO) 

and compared with the assumed value for W^. If the ratio of the new value of 
Wj to the assumed value differs from 1.0 by more than 0.001 the entire itera- 
tion procedure is repeated using the new value of W-j as the assumed value. 

This iteration is required because (AH)^ is independent of W 1 . When the ratio 
of the new value of Wj to the assumed value differs by 0.001 or less, the 

free-stream capture mass flow and capture area (both excluding inlet spillage), 
and the free-stream momentum of the capture mass flow W-j are computed by 


W o= W 1 
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In the present forebody computations the engine capture flow is penalized 
relative to friction and heat-transfer losses at the free-stream pressure. 
Physically these losses would occur at the forebody pressure and would result 
in a different entropy gain. A check into the effect of imposing the forebody 
friction and heat losses at the forebody pressure revealed an Increase in pre- 
dicted mass capture from 1.5 percent to 7 percent and an increase In thrust 
coefficient of 1.9 percent to 10 percent for the range of forebody angles 
considered. 


Inlet Flow 


T he change ip flow properties from the forebody or (or inlet entrance) 
station 1 to the 'nlet throat or station 2 is computed in a similar manner to 
that of the forebody flew. A form of the one-dimensional momentum and energy 
equations is used which accounts for Inlet friction and heat losses as well as 
inlet mass flow sp-'llage and a similar method Is utilized to approximately 
account for inlet flow shock losses. The area ratio A-j^ from inlet 
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entrance to Inlet throat, the ratio of the Inlet Invlscld kinetic energy effi- 
ciency (n K )j (assumed to be equal to .994) through the Inlet shocks, the ratio 

of inlet mass flow spillage to free-stream capture flow excluding spillage 
(W 0 " W, ) 

g and the Inlet flow heat loss (AH t ) are specified. The inlet 

o 1 

kinetic energy efficiency Is given by 
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which differs from the normal kinetic energy efficiency definition in that 
only the Inviscid inlet shock losses are included. The friction losses are 
estimated through the Reynolds analogy expression given in the following 
relationships. 
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This method of consideration of friction losses with heat losses is felt to be 
preferable to that of using adiabatic experimental kinetic energy efficiencies. 
The use of h^ Q in the Reynolds analogy expression throughout the inlet 

combustor and nozzle Is justifiable through the fact that sensitivity studies 
on varying friction losses revealed large changes in an engine's Internal 
friction losses result in small changes In engine performance. The remaining 
governing equations for inlet flow are 
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Guessed values of pressure and temperature are computed from: 


*2, S 


P 2 = P 2 
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To = T, 


As for the forebody flow computations, the thermodynamic property data of 
references 12 and 13 are used along with (p 2 » T 2 + AT) and (pg, Tg ) . The value 

of the temperature Is Improved by the equation 

S 2, S ~ S 2 
/To \ S 2 ' S 2 


T 2 " T 2 It^TaI 


where S 2 and are the entropies for the states corresponding to (pg. Tg ) 
and (pg, T 2 + AT), respectively. When successive values of Tg differ by less 

Ao 


than 1°, the area ratio 


Is computed from the relation 


p 2, S V 2, S 1 


[ 2 ( h t, 


2 " h 2, 


and compared with the value computed by equation (36). If the value of 
*2^2, S com P ute< * by equation (43) Is not within 0.00001 of the value computed 

by equation (36) the new station 2 flow values are Input Into equation (37) for 
the Isentroplc expanded flow parameters and the Iteration beginning with 
equation (37) Is repeated. 
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Combustor Flow 


1 


1 


The change In flow properties from the Inlet throat or combustor entrance 
(station 2) to the combustor exit station (station 3) Is computed utilizing a 
form of the one-dimensional momentum and energy equations. The one-dimensional 
flow equations as utilized in the combustor analysis are generated through 
application of the control volume concept on the combustor as indicated In 
figure 3. The fuel equivalence ratio 4>, combusclon efficiency n com | > » the 

combustor heat loss AH CQmb , the area ratio across the combustor A^/A^* heat 

added (J/Kg) to the combustor per kilogram of fuel ( h t)f ue i» and the tota1 heat 

added to fuel AH fue1 through the engine cooling cycle are specified. 

The one-dimensional momentum equation as obtained through application of 
the control volume concept to the combustor is given by. 


¥2 * P2 A 2 = V 3«l 




+ P3 A 3 




(pV 2 ) 


AV, comb. 
~2 


A (C f } 

H W, comb 'AV comb 


(A44) 


The wall pressure integral is given by. 



P w dA = (A 3 - a 2 ) (p 3 + p 2 ) ■£ 


(A 45) 


The parameter B p is the combustor wall pressure Integral factor of reference 

19 which relates the wall pressure integral with the value for the Integral 
assuming tne wall pressure distribution can be approximated by the arithmetic 
average of the combustor entrance (station "2") and exit (station "3") pressures. 
The form of Reynolds analogy assumed is given by. 


AH 


comb 


_ \ 


AV ; comb 


,( 


p V/A7, co mb 


t, 0 


^W, comb 


]V 


comb 


2Pr v 


(A46) 


(pV 2 ) 


A V, co mb . 

T " w, comb. 


(C c ) 

‘ AV comb 


2 

Pr^ AH comb Vv, comb 
^t, 0 h W, comb^ 


(A47) 
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t 
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Substituting equation (47) into equation (44) and rearranging gives 


*'" 2 " 


**’*-*■- ■ Ts.ry j 

w, (1 + 


(A48) 


The one-dimensional energy equation as obtained through application of 
the control volume concept on the combustor is given by 



' 4 

h 2 * T 


AH 


comb 


AH 


W, 


fuel 


<h t ) 


fuel 



' n comb' * h 'fuel 




1 


(A49) 


The use of the combustion efficiency in the manner shown in equation (49) 
accounts for the reduction in heat release corresponding to the portion of the 
fuel not burned. The one-dimensional energy relation for the combustor exit 
(station 3) gives 


y2 

h 3 * h t, 3 " T (A50) 

The method of solution for the flow conditions at the combustor exit (station 3) 
consists of making a first guess for the pressure p^ and temperature T^» 

computing a first guess for the velocity from equation (48) and computing 

a first guess for the enthalpy h, from equations (49) and (50). Thermodynamic 

property dati are then used to obtain first guess values for the remaining 
combustor exit conditions. The first guess for the combustor exit mass flow 
is then computed from 


( W 3^G ~ p 3 V 3 A c 


(A51 ) 
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and compared with the known cembdstor «ft mass flow given by 



The combustor exit static pressure is then corrected by 


(A52) 


M 3 

p 3 = (p 3*G 


(A53) 


This iteration procedure is repeated until the ratio of the computed (or guessed} 
mass flow from equation (51} to the known mass flow of equation (52) differs 
from 1.0 by 0.0001. 


Nozzle Flow 


The change in flow properties from combustor exit (or nozzle entrance) 
station 3 to the nozzle exit station 4 is computed using a form of the one- 
dimensional momentum and energy equations assuming equilibrium chemistry. 

The form of these equations also accounts for nozzle flow friction and heat 
losses. The area ratio A^/A^ from the nozzle exit to the inlet entrance and 

the nozzle heat losses are specified. 


For a well designed nozzle (no strong shock waves) the flow is essentially 
isentropic in nature and therefore the method of computation consists of first 
expanding the flow from the combustor exit (station 3) to the nozzle exit 
(station 4) isentropicaliy corresponding to the area ratio given by 


(A54) 


The total entAilny at station 4 is given b> 



h 


t, 


t, 3 


AH, 



(A55) 
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The method of solution and the Ismtropfe relations are as follows. 

Isentroplc Relations 



F-1+^M 2 


E , K-1 
L 1TTT 


Guessed values of pressure and temperature are computed from 



(A5S) 

(A57) 

(A58) 


(A59) 


(A60) 


As was for the forebody flow computations, the thermodynamic property data are 
then used along with (p 4 A , T 4 A + AT), and (p^ A , T^ A ) to obtain values 

of sj and S^. The guessed value of the temperature Is Improved by 


T 


4, A = T 4, A 



(A61) 


where A and S* A are the entropies for the states corresponding to 
( p 4 A’ T 4 and ^ p 4 A’ + respectively. When successive values of 
a differ by less than 1 ° the area ratio A^/A^ Is computed from the 
relation 
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and compared with the value computed by equation (54). If the value of A^/Aj 

computed by equation (62) Is not within 0.00001 of the value computed by equa- 
tion (54), the new flow values corresponding to the subscript (4, A) are Input 
into equation (56) for the station 3 parameters and the iteration beginning 
with equation (56) is repeated. 


The momentum relationship assumed for the one-dimensional nozzle flow is 
given by 



2 

Pr^ AH W V A V. W 
(h t, o N ) W 3 


(A63) 


and the one-dimensional energy equation assumed is given by 


V 2 v 2 AH 

h - h + 3 v 4 ah N 

h 4 " h 3 + T " 1 WT 


(A64) 


The first guess of the nozzle exit static pressure and temperature is made 
using the relations of adiabatic flow of a perfect gas in a constant area chan- 
nel with friction. The Mach number is computed by iterating for the Mach num- 
ber using the relation 


M 


a 



X A 



(A65) 


where. 


F = 1 



(A66) 


and is assumed equal to K A The value cf ^ Is used as the first 
guess for (in on right-hand side of equation (65))and the iteration is 
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repeated until successive values of M 4 differ by less than 0.001. When this 

tolerance Is satisfied the first guess static pressure and tanperature are 
computed using 



(A67) 

(A68) 


With a guess for p 4 (equation (67)) and the value for h 4 computed by 

equation (64), an iteration on the temperature is performed using as a first quess 
the value of T 4 computed by equation (68). Succeeding guesses for T 4 are 

computed from 


T 4 = (T 4 - AT) + 2AT — j (A69) 

Thermodynamic property data of references 12 and 13 are used to find the 
values of h for (T 4 - AT) and h + for (T 4 + AT), respectively, and the Itera- 
tion is repeated until succeeding guesses for T. differ by less than 1°. The 
area A 4 is computed from the relation * 



(A70) 


The area A 4 as computed by equation (70) is compared with the value of 
A 4 as computed from 



(A71) 


If the nozzle exit area computed by equation (70) is not within 0.01 per- 
cent of the value computed by equation (71) a new guess for the value of p. 
is made using 4 
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JL 

P 4 V 4 




(A72) 


The iteration is then repeated beginning^ th $uation J 69) ° nt11 th# 
values of the nozzle exit area A 4 a*ree within 0.01 percent. 
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TABLE I 


Parameters Needed for Engine 
Performance Computations 


Free Stream 


1. Altitude (meters) 

2. Free-stream Mach number, M 

3. (h t, o ' h W ) ( Joules / K 9) 

4. Guessed captur° mass flow (Kg/sec) 

A c - A 1 

5. Spillage flow parameter 

M c 


Forebody Flow 

1. Pressure ratio across forebody shock, P-|/P 0 

2. Inviscid forebody flow Mach number, M^ Iny 

3. Kinetic energy efficiency across forebody shock, [\)f 

4. Loss of heat, by capture flow to forebody, (AH) f (Ooules/sec) 

5. Cowl area T 

6 . Forebody angle to flight direction, w (deg) 


Inlet Flow 


1. Inlet contraction ratio, A.j/A 2 

2. Kinetic energy efficiency across inlet shocks, (n^)j 

3. Loss of heat by capture flow to inlet surfaces, (AH)j (Joules/sec) 


Combustor Flow 


1. r uel Equivalence ratio, <t> 

2. Combustion efficiency, (n) com jj 

3. Combustor exit area to combustor entrance area, (A 3 /A = 0.64) 

4. Combustor wall pressure integral parameter. Bp 

5. Loss of heat by capture flow to combustor surfaces, (AH) com j ) (Joules/ sec) 

6 . Heat added to the combustor per kilogram of fuel, (^fuel (J°ules/Kg) 

7. Heat adced to fuel due to regenerative cooling of engine, (AH)-,, , 

(Joules/ sec) Tue 
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TABLE I.- Continued 


Nozzle Flow 

1. Equilibrium flow is assumed 

2. Area ratio across nozzle, A^/A^ 

3. Nozzle kinetic energy efficiency, (n^ 

4. Loss of heat by capture flow to nozzle surfaces, (AH) N (Joules/sec) 

5. Nozzle exit flow angle with forebody surface, B (deg) 

Stream Tube Force Parameters 

1. Pressure force on the capture flow stream tube over the spillage area 
and normal to forebody surface, Aj (N) 

2. Pressure force on the lower surface of the plume normal to forebody 

surface, A (N) 

3. Inlet additive drag, (N) 

4. Nozzle olume drag, (N) 
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Figure 1.- Inner module of sc ramjet engine concept; H = 0.46 meters and 
sweep angle equals 48 degrees in the present engine. 
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(a) Pressure ratio across vehicle forebody shock. 
Figure 4.- Nonviscous forebody flow. 
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(b) Fngine module capture flow. 
Figure 5.- Concluded. 





Figure 6.- Inlet aerodynamic contraction ratio (experimental). 
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Relationship of present engine with vehicle assumed (for engine 
performance computations) including nozzle-vehicle afterbody 
relationship (dimensions in meters). Cowl external surface 
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Figure 9.- Plume drag coefficient. 
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e 10.- Plume force coefficient 










Figure 11.- Nozzle exit area. 
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Figure 14.- Capture flow heat losses for M = 6.0. 





it transfer ratioed to heat transfer 
inferred from reference 5. 
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Figure 16.- Capture flow heat losses. 
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Figure 16.- Concluded. 
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iheoretical nozzle exit flow area for the heat sink engine 
of reference 18 ratioed to the nozzle exit flow area for 
<t> = 1.0 as compared with the same parameter of the 
present engine. 
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Figure 19.- Theoretical mass weighted nozzle exit flow anqle for 
the heat sink engine of reference 18 ratioed to the 
nozzle exit flow angle for 0 = 1.0 as compared with 
the same parameter for the present engine. 





Theoretical plume drag coefficient for heat sink engine 
of reference 18 ratioed to the plume drag coefficient for 
<D = l.O as compared with the same parameter for the 
present engine. 
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i the flight direction. 

21.- Concluded. 











re 22. - Thrust coefficient of Langley scramjet engine. 






Figure 23. - Specific impulse of Langley scramjet engine. 






Effect of vehicle forebody angle on mass capture and force 
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Figurfi 25.* Effects of engine capture mass flow heat losses with engine friction 
losses held constant. 
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Figure 28.- Additive drag, D. , and nozzle plume drag. D p , ratioed to engine 
force A. 
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